Prom the Fermi glass towards the Mott insulator in one dimension: 
Derealization and strongly enhanced persistent currents 
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When a system of spinless fermions in a disordered mesoscopic ring becomes instable between 
the inhomogeneous configuration driven by the random potential (Anderson insulator) and the 
homogeneous one driven by repulsive interactions (Mott insulator), the persistent current can be 
enhanced by orders of magnitude. This is illustrated by a study of the change of the ground state 
energy under twisted boundary conditions using the density matrix renormalization group algorithm. 
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For a fixed kinetic energy (t = 1), the interplay be- 
tween the disorder (W) and a repulsive interaction (U) 
is one of the central problems of mesoscopic physics [Q . 
For instance, its understanding appears as a necessary 
step towards explaining why the measured persistent cur- 
rents Q J in mesoscopic rings are about two orders of 
magnitude larger than the value predicted by the non- 
interacting theory jjj. This discrepancy has prompted 
to study the role of the interaction using perturbation 
theory assuming a Luttinger liquid and using renor- 
malization group |H for the disorder, or considering the 
limit where the rotational invariance of the interacting 
system is broken by the disorder || . On the other hand, 
exact numerical calculations on small Id models of spin- 
less fermions have led to the conclusion that repulsive 
interactions cannot significantly M enhance J, being on 
the contrary rather detrimental [g| to the expected ef- 
fect. It was then argued that this is an artifact of spinless 
Id models, and that the full problem with spin does ex- 
hibit H the expected enhancement. We show in this Let- 
ter that even without spin, repulsive interactions can give 
rise to a large enhancement associated with a charge reor- 
ganization of the ground state. The interaction for which 
the ground state is reorganized fluctuates from sample to 
sample. Thus, the effect could not be seen in the previ- 
ous studies of the ensemble averages, but it becomes very 
striking if one studies the orbital response of individual 
mesoscopic samples as a function of U. 

Another motivation to re-examine the ground state of 
one dimensional spinless fermions stems from Shepelyan- 
sky's proposal concerning the role of interactions in an 
Anderson insulator J^]: two interacting particles (TIP) 
can be localized on a length L2 much larger than the 
one-particle localization length L\. Large ratios L2/L1 
characterize the states in the bulk of the TIP-spectrum. 
However, a TIP ground state built from two one-particle 
states localized far away from each other cannot be reor- 
ganized by a short range interaction. Interesting effects 
on the ground state properties require reasonably high 
filling factors. For the finite density system, a Fermi liq- 



uid approach was suggested by Imry (!(]] and later devel- 
oped in Refs. (ll]jlj], reducing the problem to the study 
of the derealization of a few quasi-particles above the 
Fermi sea. From a scaling argument, one obtains that 
two quasi-particles need an excitation energy e«( (irrel- 
evant for the usual applications in solid state physics) in 
order to recover the enhancement factor L2/L1 of two 
(excited) bare particles. This quasi-particle approach 
may be appropriate if there is no dramatic reorganiza- 
tion of the ground state by the interaction. We want to 
consider another limit: when the Fermi liquid becomes 
instable as the Mott insulator is approached. To avoid 
any confusion, let us say that the Fermi liquid will be bet- 
ter named a Fermi glass in the disordered limit which we 
consider, characterized by Anderson localization without 
interaction. By Mott insulator, we mean a strongly cor- 
related array of charges which is pinned and distorted by 
the random lattice. 

For the TIP problem (symmetric states of two elec- 
trons with opposite spins and on-site interaction), the 
ratio L 2 /Li is maximum fl3| when the system becomes 
instable between two limits: free bosons for U = and 
hard core bosons for U — oo. When U « t, the mix- 
ing of the one-particle states by the interaction is maxi- 
mum and the system becomes weakly chaotic with criti- 
cal spectral statistics for system sizes Mwli. This crit- 
ical regime cannot be described by perturbation expan- 
sions in powers of U or t 2 /U, starting from the free boson 
and hard-core boson limits respectively. Therefore, ap- 
proaches starting from simple limits for the ground state 
(Luttinger liquid for weak disorder Mott insulator 
for strong interaction, Fermi glass for weak interaction 
and large disorder) may not be valid for U~t and large 
disorder, say k{l f=a 1 (k^ 1 characterizes the spacing be- 
tween the electrons in the clean system and in Id the 
elastic mean free path is l = L±). 

This suggests us to study the many body ground state 
in the most chaotic case, between the limits which corre- 
spond to the free and hard-core particles in the archetyp- 
ical TIP example. If M denotes the number of sites of 
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a Id lattice with N spinless fermions, we expect a maxi- 
mum delocalization when the parameters (U, t, W 1 N and 
M) are chosen in such a way that the system becomes in- 
stable between two insulating limits with different charge 
configurations: Anderson insulator in the Fermi limit and 
Mott insulator in the strongly correlated limit. The An- 
derson insulator yields an inhomogeneous charge config- 
uration: strongly localized states are populated with the 
restriction imposed by Pauli principle. The Mott insula- 
tor yields a more or less homogeneous charge configura- 
tion. The crossover between those two limits will yield 
a profound spatial reorganization of the ground state, 
making the system more sensitive to any external per- 
turbation as a twist in the boundary conditions or a flux 
$ in a ring. 

To illustrate this, we consider spinless fermions on a 
disordered chain with nearest neighbor (NN) interaction 

M M M 

H = -ty"](c|cj_i + c|_ 1 Ci) + Villi + n i n i-l 
i—1 i—1 i—1 

and twisted boundary conditions, cq = exp(i$)cA/- The 
operators Cj (cj ) destroy (create) a particle on site i and 
rii = c|cj is the occupation operator. The on-site ener- 
gies Vi are drawn from a box distribution of width W . 
Setting t = 1, the strength of the disorder W and the 
interaction U are measured in units of the kinetic en- 
ergy scale. Though strong long-range interaction always 
imposes an homogeneous configuration, we restrict the 
study to interaction between NN. This is why we begin 
to discuss half filling, where the ground state will be a pe- 
riodic array of charges located on the even or odd sites of 
the chain when U — > oo. This is the homogeneous config- 
uration driven by the interaction. On the contrary, when 
U — > 0, an inhomogeneous configuration is obtained for 
a sufficiently strong disorder. 

The numerical results are obtained with the density 
matrix renormalization group (DMRG) algorithm [fL5| . 
The calculation of ground state properties of interacting 
fermions in disordered one-dimensional systems is 
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FIG. 1. Charge configuration for a typical sample (d of 
Fig. |) for TV = 10 particles on M = 20 sites at W = 9. 
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FIG. 2. Phase sensitivity D(U) for four different samples 
with N = 10, M = 20 and W = 9 in (decimal) logarithmic 
scale. W and U are measured in units of t. Thick dots and 
inset: average of log(D). 

possible with an accuracy comparable to exact diagonal- 
ization, but for much larger systems (where we keep up 
to 2000 states per block fig]). 

The reorganization of the ground state induced by the 
NN repulsion is shown in Fig. [lj where the density p (ex- 
pectation value of rii) is plotted as a function of U and 
site index i. To favor the inhomogeneous configuration, 
the disorder is taken large (W = 9) and L\ 100/M^ 2 is of 
order of the mean spacing fcjT 1 = 2 between the charges. 
For U «s 0, one can see a strongly inhomogeneous and 
sample dependent density, while for large U a periodic ar- 
ray of charges sets in. These two limits are separated by a 
sample dependent crossover regime. For certain random 
configurations, the periodic array is quickly obtained by 
a weak repulsive interaction, while one needs a strong 
interaction for other samples. 

To measure the delocalization effect associated to this 
change of configuration, we study the phase sensitivity of 
the ground state. The energy difference between periodic 
($ = 0) and anti-periodic ($ = n) boundary conditions, 
AE = (— ) N (E(0)— E(ir)) conveys similar information, in 
the localized regime, as other measures of the response 
of the ground state to an infinitesimal flux threading the 
ring: the Kohn curvature (charge stiffness) oc E (<I> = 0) 
and the persistent current J cx — E ($ = 0). For strictly 
lei systems, the sign of E(0) — E(tt) simply depends on 
the parity of N, and the factor (— ) N makes AE positive. 

The phase sensitivity, defined by D(U) = (M/2)AE, 
is shown in Fig. || for four samples at half filling with 
W = 9. Both for C7«0 and [7>1, D(U) is exponentially 
small, but sharp peaks appear at sample dependent val- 
ues U c , where D(U C ) in certain samples can be 4 orders 
of magnitude larger than for free fermions. Remarkably, 
the curves for each sample do not present any singularity 
at U — which could allow to locate the free fermion case. 
Peaks can be seen at different sample dependent values 
of U (positive or negative) . Attractive interactions favor 
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FIG. 3. Phase sensitivity for the sample d of Fig. |^. Solid: 
D(W) at 17 = 0.9; dashed: D(U) at W = 2. Inset: Probability 
distribution (dots) of rj = log D([7 = 2) - log D((7 = 0) calcu- 
lated from 10000 samples (M = 20, N = 10, W = 9) fitted by a 
Gaussian with variance <r 2 = 1.46. 

the more inhomogeneous density. Repulsive interactions 
favor the more homogeneous density. Free fermions cor- 
respond to an intermediate case. For small repulsive in- 
teractions, the system is an Anderson insulator slightly 
delocalized by U, and D(U) increases as a function of 
U. At U c , the regular array of charges is established, 
and thereafter it becomes more and more rigid (pinned 
by the random lattice); thus D(U) decreases as a func- 
tion of U. Certain samples go from the inhomogeneous 
density to the periodic array in a few steps signaled by 
additional peaks of sensitivity. Examining the U depen- 
dence of the density of those samples, one can note local 
defects in the periodic array subsisting up to large values 
of U. The thresholds U c are strongly sample dependent 
giving rise to a very wide distribution of phase sensi- 
tivities: the ensemble average at a given U mixes very 
different behaviors and provides rather misleading infor- 
mation. As shown in Fig. |^, (log D(U)) decreases for 
repulsive interactions, except for a small interval around 
U ~t (inset) where a local maximum is obtained. Using 
5000 samples we confirm beyond statistical uncertainty 
that repulsive interactions delocalize in certain parame- 
ter region. However, this small average effect is not rep- 
resentative of the dramatic enhancement characterizing 
individual samples. 

We obtain approximately log-normal distributions for 
D(U) as well as for the parameter r\ — log D(U = 2) — 
log D(U = 0) that measures the relative increase of the 
charge sensitivity with respect to the free fermion case. 
The width of the 77-distribution depends on U, and for 
U = 2 we can see that variations of D over more than an 
order of magnitude are typical (inset of Fig. ||). 

For weak disorder (W — 2,Li M) we recover 
the expected behavior starting from the clean limit, 
using bosonization and renormalization group argu- 
ments Jl6|,[f4| : a repulsive interaction reinforces localiza- 
tion, in contrast to a (not too strong) attractive interac- 
tion which delocalizes (Fig. ||). Fixing U = 0.9, we show 
(in the same graph) how one goes from the weak to the 
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FIG. 4. Phase sensitivity for M = 40 and N = 10 for 3 
samples with W = 4 (thin lines) and W~5 (thick lines). 

strong disorder limit in a given sample (the same seed of 
the random number generator is kept). The phase sen- 
sitivity decreases when we increase the strength of the 
potential fluctuations, except around W — 9 (kfl « 1) 
where the charge reorganization takes place. The con- 
clusion that a repulsive interaction favors localization is 
no longer valid. We underline that we have considered 
in this study larger disorder than in Ref. |l6). The ef- 
fect is reminiscent of the Coulomb blockade phenomenon. 
When the occupation numbers are approximately good 
quantum numbers (0 or 1) transport is blocked; in the 
transition between such extreme configurations the oc- 
cupation numbers are no longer good quantum numbers 
and transport is favored. 

Now, we study if a repulsive interaction can delocalize 
at a weaker disorder for a lower filling. Of course, if wc 
keep NN interactions, U will no longer yield a fully pe- 
riodic array (in contrast to a longer range interaction): 
the ground state can partly match the random potential 
without having NN electrons. One expects a more local 
and partial ground state reorganization than the global 
one discussed for half filling. In Fig. |], we show some 
typical behavior obtained at quarter filling for W — 4 
and 5. For certain samples (dash-dotted), the U = 
Fermi glass does not have NN occupied, the charge con- 
figuration does not change when turning on the interac- 
tions, and D(U) decreases with U. For others we can 
observe local charge reorganization accompanied by an 
increase of D(U) at one (dashed) or two (solid) values 
of U. Higher disorder within the same random configu- 
ration (thicker lines) necessitates a stronger interaction 
U c to produce the charge reorganization. In our model 
with NN interactions and high disorder, the increase of 
the phase sensitivity is related to the occurrence of NN 
in the Fermi glass. Elementary combinatorics dictates 
that for a given filling factor x — N/M and large M, the 
probability of obtaining a configuration without NN is 

p „ e Mg{x) ( ) = j ( C 1 - x) 2{1 ~ x) \ 

Since g(x) <0 in the interval (0, 1/2) of interest, configu- 
rations with NN occur with probability 1 in the large M 
limit. 
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FIG. 5. Top: phase sensitivity of three samples with 
W = 2, N = 8 and various fillings: M = 20 (solid), M = 40 
(dotted), and M — 80 (dashed). Bottom: charge density at 
U = (solid) and (7 = 10 (dashed) for the sample with A/ = 80 
that shows a peak around U = 1.5. 



(Mott). When fcfi < 1, the particles can be strongly lo- 
calized far away from each other, and a short range inter- 
action does not affect a strongly insulating ground state 
(Fermi glass). Only very excited quasi-particles can be 
delocalized by the interaction. Between those two limits, 
kfl re 1, the ground state can be deeply reorganized by 
a repulsive interaction [/ret. This reorganization is ac- 
companied by a large derealization effect. Preliminary 
results Jl7| for the interaction induced charge reorganiza- 
tion in two dimensions lead to similar conclusions. This 
might be related to the problem of the metal-insulator 
transition jl8j observed in Si-Mosfets at very low fillings. 
In quasi-one dimension, this might give some insight for 
the persistent current of a disordered multichannel ring. 

We thank G. Ingold, J.-F. Joanny, I. Safi, X. Waintal 
and H. Weidenmiiller for useful discussions. Financial 
support from the TMR network " Phase coherent dy- 
namics of hybrid nanostructures" of the EU is gratefully 
acknowledged. 



Fig. H presents D(U, x) at a smaller disorder (W = 2) 
for filling factors decreasing down to 1/10 (N = 8 and 
M = 20, 40, 80). We underline that such a study is clearly 
out of reach if one is restricted to standard diagonaliza- 
tion algorithms. To estimate a sample dependent length 
L x , one assumes that D(U = 0) re (M/2) exp(-2M/L x ). 
This gives for M = 20(40) L\ re 17(10), corresponding to 
kfl > 1. In this case a repulsive interaction localizes. On 
the other hand, for a larger size M = 80, one recovers the 
derealization effect when D(U = 0) re 10~ 6 , i.e. when 
L\ re9re M/N (kfl re 1). The peaks can be now broader 
and the magnitude of the effect is weaker (two orders of 
magnitude). Studying the charge density (bottom), we 
notice a slow local charge reorganization induced by a 
large variation of U. When D(U) re 10 -8 (one of the 
samples with M = 80), the scale Li re 7 of those sam- 
ples is smaller than M/N, kfl becomes too small and the 
Fermi glass cannot be reorganized by a NN interaction 
at this low filling. 

We summarize the main conclusions that we draw from 
this DMRG study of one dimensional spinless fermions. 
(i) Each sample should be individually studied. The 
(log)-averages over the ensemble are not representa- 
tive, (ii) The ratio L\/{M/N) (say kfl) defines differ- 
ent regimes for short-range repulsive interactions. For 
kfl > 1, the interaction establishes a correlated array of 
charges inside L\ which is pinned by the random lattice. 
The larger is U , the more rigid is the array, the more effi- 
cient is the pinning and the system is strongly insulating 
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